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THE CALCULATION OF N- DIMENSIONAL 
PROBABILITY ELLIPSOIDS 

E. R. Lancaster 
Theoretical Mschanics Branch 

ABSTRACT 

The problem is to find a real number c such th»?it the probability 
is p that a normally distributed random variable x will lie in an 
ellipsoid defined x^Ax = c, where A is the inverse of the covariance 
matrix of x. This problem arises, for example, in orbit determination 
when it is desired to know the dimensions of an ellipsoid centered about 
+he nominal value of the state vector at a given time such that the 
probability is p that the state error vector will be in the ellipsoid. 


THE CALCUUTION OF N- DIMENSIONAL 
PROBABILITY ELLIPSOIDS 


INTRODUC!HON 

If y is a real normally distributed random vector, the probability 
that y is in a region S is 


where 


p » f(x)do) 
f(x) - ft e-^ A * 


( 1 ) 

( 2 ) 


and 


du) = dXj^dxa . . .dx^, (5) 

ft =. |A|^/(2rr)^. (4) 


The real vectors x and y have n components, x* is the transpose 
of X, and A is an n-by-n symmetric positive-definite matrix of real 
numbers with determinant |a| . 

The equation 

x'Axwc, ( 5 ) 

where c is real, defines an ellipsoid in n- space with center at 
X ss 0. The problem is to find c, given that p is hnovm and that S 
is defined by ( 5 ) • The matrix A is the inverse of the co-variance 
matrix of y. The mean of y has been talten to be 0. 

This problem arises, for example, in orbit determination when it 
is desired to know the dimensions of an ellipsoid centered about the 


1 



nominal value of the state vector at a given tine such that the 
prohahillty is p that the state error vector will he in the ellip- 
soid. Instead of the full state error vector y we can have the error 
vector for a suhspace of the entire state space, i.e,, a vector of 
dimension less than n whose components form a subset of the full 
error vector. Cils is true since the marginal distribution of any set 
of components of y is normal with mean 0 and variances and co- 
variances obtained by taking the proper con^onents of the inverse of A. 

INTEGRALS OF QUADRATIC FORMS 

As a step toward solving the problem described in the introduction, 
we consider integrals of the form 

g(A) =* f f(x' A x)du), (6) 

where A and do) are defined in the introduction and S is the ellipsoid 
of (5). Since A is positive definite, there exists a nonsingular 
matrix C such that 

A c » I, (7) 

where I is the identity matrix and C* is the transpose of C. 

Let 

X =? C z . (8) 

Then x’ Ax»z’ C’ AC z = z*z. 

The Jacobian of the transformation (8) is C, the absolute value of the 
determinant of C. Thus (6) becomes 
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(9) 


g(A) - C J f(z*z)dc 

where da « dz^ dZ 2 . . . dz^, and R is the sphere with center at 
defined hy 

z*z ■ c. 


We now transform to generalized spherical coordinates by 


Zj *» r cos 9 ^ 
i-1 

”(ir sin cp.^ cos 9^; i 


J=1 

n-2 


2; 5> . . . ^ n-2 


' (T[ sin 9 j^ cos 0 


z = r 
n 


n-a 

(tt 

>1 


sin 9j^ sin 0 


where 0 sC 9 ^CTf; 0^r<®,n^5* n«2, z^ » r cos 0, 

zg a r sin 0 . The Jacobian of this transformation is 


n-2 


J = r“"^ I I 9,, a a J. 

J =1 ^ 


Also 


z*z = r^. 


If n»2, J^r • 

Thus (8) becomes 


c 2TT 

^ .n-1 


g(A) =» C r- f(r^)dr d0 
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» 0 

do) 


Prom (7) 


c - i/|a|* , 


n-2 TT 

J tp d«p - r"'^ (i)/r(n/2), 


J=1 'o 


where r is the gamma function, ^us we have 


g(A) - C K f. r"’^ f(r®)dr, 

W 


( 11 ) 


where 


K-2ITr“"^(i)/r(n/2), 


l.e., K = srr* /(*-!)! If n = 2A, 

K = 2 **^ tl*/[l*3*5'7”*(a*-l)] if n = 2 A + 1 , 

where ft is an integer s 1. 

THE NORMAL PROBABILITY ELLIPSOID 
Applying (ll) to (l) gives 
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6r, n » 2-ft, 


(12) 


"3 = 


2//ff r 

l-5*5---(2ft-l) J, 


2ft -i r^a 
r e ® dr, 


n » 2-ft + 1. 


(13) 


These integrals can he evaluated hy formula 313 -5 of reference 1. 
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The function p(c) Is monotone increasing on the interval (O^ ») and has 
a single Inflection point in the interval. Ihus the Kievton-Bapheon 
method will find the value of c for a given value of p If we use the 
value of c correspcoiding to the inflection point as the starting value. 
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